  x = [0  .999999999999844  -.499999558846537  -.499999906019974  -1.35808453673367e-007  ];

  y = [0  0.  -.866025658484298  .866025458043624  5.29764504205778e-007  ];


  z = [1  5.58120440369433e-007  3.257232895128e-007  -5.91106152408667e-007  -.99999999999985  ];
  
 % plot3(x,y,z,'rs','Markersize',15);
 hold on;
 % sphere;
 % axis equal;
 % hold on;
  plot3(x,y,z,'rs','Markersize',10,'Markerfacecolor','blue');
  hold off;
  
 energy = 0;
for i = 1:length(x)
    for j = i+1:length(x)
        s = sqrt((x(i)-x(j))^2+(y(i)-y(j))^2+(z(i)-z(j))^2);
        energy = energy + 1/s - 1/2*(log(s) + log(2+s));
    end
end